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For  the  low  temperature  applications  of  Bi2Te3  based  thermoelectric  power  generators  (TEGs),  determi¬ 
nation  of  low  temperature  material  behaviors  is  important.  In  the  temperature  range  of  100-375  K,  tem¬ 
perature  dependency  of  Seebeck  coefficient  and  electrical  conductivity  of  a  Bi2Te3  based  TEG  are 
experimentally  determined.  Furthermore,  for  a  constant  temperature  difference,  the  variation  of  maxi¬ 
mum  power  output  with  the  mean  temperature  is  analytically  examined  based  on  the  experimentally 
measured  material  properties.  It  is  observed  that  250  K  seems  a  critical  mean  operating  temperature 
for  the  considered  Bi2Te3  module.  Correlations  for  temperature  dependency  of  material  quantities  are 
also  given  for  the  temperature  region  of  100-375  K.  The  results  can  be  used  for  the  low  temperature 
applications  of  Bi2Te3  based  TEGs. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  devices  consist  of  simple  thermocouples  which 
can  directly  convert  either  thermal  energy  into  electrical  energy 
for  power  generation  applications  or  electrical  energy  into  thermal 
energy  for  heating  or  cooling  applications.  Although  thermoelectric 
devices  have  many  advantages;  they  are  not  commonly  used  due 
to  their  lower  efficiency  and  higher  costs  in  comparison  with  those 
of  the  conventional  systems.  After  the  recent  advances  in  material 
science,  however,  thermoelectric  power  generator  modules  with 
5%  efficiency  and  Peltier  coolers  with  COP  value  of  1.5  for  20  °C 
temperature  difference  are  commercially  available.  On  the  other 
hand,  higher  efficient  (up  to  17%  for  200  °C  temperature  difference) 
TEG  modules  based  on  different  type  of  nano-materials  have  been 
developed  in  many  research  laboratories  [1-3,19]. 

In  most  of  the  applications,  thermoelectric  devices  are  used  for 
cooling  purposes  [4-7].  However,  thermoelectric  generators  have 
also  been  used  in  a  number  of  applications  such  as  battery  charg¬ 
ing  [8],  waste  heat  recovery  [9,10],  power  from  radioisotopes 
[11,12],  as  well  as  electrification  of  rural  homes  [13].  Thermoelec¬ 
tric  devices  are  commonly  produced  for  the  operating  tempera¬ 
tures  ranging  from  room  temperature  to  higher  temperatures. 
Therefore,  there  are  only  a  few  works  in  literature  on  the  charac¬ 
teristics  and  performance  of  thermoelectric  devices  at  low  temper¬ 
ature  range  [14-17].  Heat  energy  transfer  from  environment  to  the 
liquefied  natural  gas  (LNG)  during  its  evaporation  can  be  used  to 
generate  electricity  by  TEGs  [14].  Therefore  evaporation  process 
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of  liquefied  gases  constitutes  a  possibility  for  the  low  temperature 
applications  of  TEGs. 

In  this  work,  a  Bi2Te3  based  thermoelectric  power  generator  is 
considered  since  they  are  commercially  available  and  relatively 
inexpensive.  Its  material  characteristics  such  as  Seebeck  coefficient 
and  electrical  conductivity  are  experimentally  determined  for  low 
temperature  region  (100-300  K).  Furthermore  the  variation  of 
maximum  electrical  power  output  for  a  given  temperature  differ¬ 
ence  with  the  mean  operating  temperature  is  analyzed  by  using 
the  experimentally  measured  material  properties  in  some  analyti¬ 
cal  expressions.  The  results  represent  the  low  temperature  charac¬ 
teristics  of  Bi2Te3  based  thermoelectric  materials.  They  can  be 
useful  for  the  studies  on  low  temperature  engineering  applications 
of  TEG  modules  as  well  as  the  material  science. 

2.  Experimental  set-up 

The  schematic  diagram  of  the  experimental  set-up  is  shown  in 
Fig.  1.  The  set-up  consists  of  thermoelectric  module  which  is 
clamped  between  two  copper  plates  contacting  hot  and  cold  sur¬ 
faces.  The  thermoelectric  module  is  surrounded  by  XPS  insulator 
to  reduce  heat  losses  to  ambient.  For  the  measurements,  tempera¬ 
ture,  heat  flux,  current  and  voltage  sensors  are  used.  Uncertainties 
of  voltmeter,  ampere  meter,  thermometer  and  heat  flux  meter  are 
±0.04%,  ±1%,  ±0.5%  and  ±3%  respectively.  K-type  thermocouples  are 
used  for  the  temperature  measurement  of  both  hot  and  cold  side 
lead  surfaces  of  TEG  module  as  well  as  the  interfaces  between 
the  copper  plates  and  hot  and  cold  plates.  The  heat  flux  meters 
are  settled  between  copper  plates  and  TEG  module  surfaces.  The 
system  also  includes  the  liquid  nitrogen  vessel  at  77  K  and  electri¬ 
cal  heaters  to  adjust  the  cold  and  hot  surface  temperatures,  TL  and 
Th  respectively. 
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Fig.  2.  Schematic  cross-sectional  view  of  TEG  module. 


A  commercially  available  Bi2Te3  based  thermoelectric  generator 
is  chosen.  It  generates  14.1  W  electrical  power  for  hot  and  cold  side 
temperatures  of  250  °C  and  50  °C  respectively.  The  overall  sizes  are 
56  mm  width,  56  mm  length  and  4.3  mm  height  (or  thickness).  The 
number  of  p-n  semiconductor  couples  is  Npn  =  241.  For  both  p-type 
and  n-type  single  semiconductors,  base  sizes  are  1.7  mm  x  1.7  mm 
and  its  thickness  is  1.1  mm.  The  cross-sectional  view  as  well  as 
thicknesses  of  materials  are  seen  in  Fig.  2.  Open  circuit  voltage 
and  short  circuit  current  are  8.4  V  and  6.7  A  respectively  for  TL 
=  50  °C  and  TH  =  250  °C. 


3.  Characterization 

During  the  experiments,  the  temperature  difference  between 
hot  and  cold  surfaces  of  TEG  module  is  kept  constant  at  20  K.  For 
AT  =  20  K  temperature  difference;  generated  voltage  V,  current  J, 
heat  fluxes  at  hot  and  cold  sides  (q^  and  q'[)  as  well  as  tempera¬ 
tures  T  are  measured  for  different  electrical  load,  RL  values  practi¬ 
cally  ranging  from  zero  to  infinite  ohms.  The  ambient  temperature 
is  set  to  25  °C  (298  K)  and  the  heat  exchange  between  ambient  and 
the  lateral  walls  of  the  generator  is  neglected.  Because  the  lateral 
surfaces  are  insulated  by  XPS  material  and  also  the  total  lateral 
surface  area  is  too  small  in  comparison  with  hot  and  cold  side  sur¬ 
faces.  Measurements  are  repeated  for  different  mean  tempera¬ 
tures,  T  =  (TH  +  TL)/2,  ranging  from  100  K  to  375  K.  Since  the 
temperature  difference  is  small  enough  (AT  =  20  I<),  temperature 
dependence  of  Seebeck  coefficient  and  electrical  conductivity  are 
assumed  to  be  linear  for  this  small  temperature  interval  ranging 
from  Tl  =  T  -  AT/2  to  TH  =  T  +  AT/2.  In  this  case,  integral  aver¬ 
ages  of  these  quantities  for  temperature  interval  of  20  K  are  equal 
to  their  values  at  mean  temperature.  Therefore,  open  circuit  volt¬ 
age  V0  is  expressed  as 


rTH  rT+AT/2 

Vo  =  Npn  /  a np(J)dT  =  Npn  /  (Aa  +  BJ)dT 

Jtl  JT-AT/2 

=  Npn  (Aa  +  BJ)AT  =  Npn0Cnp(J)AT.  (1) 

where  otnp  =  otp-  otn  is  the  Seebeck  coefficient  of  a  single  p-n  couple, 
Aa  and  Ba  are  the  linear  coefficients  representing  the  temperature 
dependency  of  the  Seebeck  coefficient  for  a  small  temperature 
interval.  Therefore  temperature  dependent  Seebeck  coefficient  is 
simply  calculated  by  measuring  the  open  circuit  voltage  and  tem¬ 
peratures  of  hot  and  cold  sides  as 

MT>-N^r  12) 

The  full  derivate  of  a  material  quantity  (Seebeck  coefficient  and 
electrical  conductivity)  with  respect  to  the  measured  quantities 
represents  the  uncertainty  of  the  material  quantity  as  a  function 
of  uncertainties  of  the  measured  quantities.  By  using  the  linear 
approximation  (small  uncertainties),  the  ratio  of  the  full  derivative 
of  the  material  quantity  to  its  value  is  assumed  to  be  the  value  of 
total  percent  uncertainty.  Therefore,  total  percent  uncertainty  of 
the  Seebeck  coefficient  is  calculated  from  (2)  as  follows 

da_dVo_dAT 

a  _  V0  AT  '  U 

and 


Aa  —  d=A  v  ±  A  at- 


(4) 


Uncertainty  of  the  measurement  of  Seebeck  coefficient  is  in  or¬ 
der  of  ±0.54%. 

Determination  of  internal  resistance  based  on  simultaneous  DC 
current  and  voltage  measurements  causes  wrong  results  since  the 
DC  current  itself  induces  a  Peltier  effect  which  changes  the  temper¬ 
ature  difference  during  the  measurement.  This  induced  tempera¬ 
ture  difference  changes  the  net  electrical  potential  difference 
applied  to  the  material  and  causes  wrong  resistance  measurement. 
Therefore  to  prevent  the  induction  of  Peltier  effect  during  the  mea¬ 
surement  of  internal  electrical  resistance,  it  is  measured  by  an 
impedance  meter,  which  uses  AC  signal.  Internal  resistance  can 
be  expressed  as 


Ri  =  Npn 


1  dn  1  dp  2  ^con 

On  An  Op  Ap  Acon 


(5) 


where  a  is  the  electrical  conductivity,  S  is  the  thickness,  A  is  the 
cross  sectional  area,  rcon  is  the  contact  resistance  between  metal 
electrodes  and  semiconductor  p-n  couple  while  the  subscripts  n 
and  p  represent  n-type  and  p-type  semiconductors  respectively. 
Here  the  resistance  of  metal  electrode  itself  is  neglected  since  the 
electrical  conductivity  of  metals  is  too  high  in  comparison  with 
those  of  semiconductors,  oe  >  {crn,  ap}. 

In  macroscopic  thermoelectric  devices,  electrical  contact  resis¬ 
tance  between  the  semiconductor  and  metal  electrodes  has  been 
reported  to  typically  be  between  10-8  and  10-9  Q  m2  [18].  In  this 
work,  rcon  is  assumed  to  be  5.10-9  Cl  m2  as  an  average  value.  Since 
Sn  =  dp  =  5pn  and  Acon  =An=Ap  =  Apn  for  the  considered  TEG,  (5)  can 
be  simplified  as 


Ri^NJL+L+^)^  =  NJ^±5z+?2m\^L 


dpn  )  A^ 


OnOD 


dpn  J  Ap] 


=  nJ—  +  ^]^. 


dpn  J  Ap j 


(6) 


Therefore  electrical  conductivity  of  p-n  pairs  is  expressed  by 


Vpn(T) 


OnOp 
On  +  Op 


Ri(T)  Apn  2 XCOn 

Npn  dpn  dpn 


(7) 
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Uncertainty  of  the  measurement  of  electrical  conductivity  is  in 
order  of  ±3.5%.  It  should  be  noted  that  the  temperature  dependen¬ 
cies  of  electrical  conductivity  are  assumed  to  be  linear  for  small 
temperature  interval  of  20  K  during  the  measurements.  Therefore, 
integral  averages  of  this  quantity  for  20  K  interval  are  equal  to  their 
values  at  mean  temperature,  like  the  Seebeck  coefficient  given  by 
(2).  To  make  the  results  more  general,  dimensionless  Seebeck  coef¬ 
ficient  and  electrical  conductivity  are  defined  as  follows 

a(T)  =  ocnp(T)/ocnp(T  =  300  K),  (8) 

a(T)  =  o nP(T) / o np(T  =  300  K).  (9) 


4.  Experimental  results  and  discussion 

During  the  experiments  ambient  temperature  is  kept  at  25  °C 
(298  K)  and  the  heat  exchange  between  ambient  and  the  lateral 
walls  of  TEG  is  neglected  due  to  the  reasons  explained  in  Section 
3.  Second  order  polynomial  expressions  are  fitted  to  the  experi¬ 
mental  data  of  the  dimensionless  Seebeck  coefficient  and  electrical 
conductivity. 

In  Fig.  3,  it  is  seen  that  the  Seebeck  coefficient  of  a  Bi2Te3  p-n 
couple  decreases  with  decreasing  temperature,  up  to  80%  for 
T=  100  K  due  to  the  lack  of  thermal  excitations  at  lower  tempera¬ 
tures.  Solid  lines  represent  the  expected  values  while  the  dots 
stand  for  experimental  values.  The  mean  absolute  percentage  error 
(MAPE)  of  the  solid  line  is  defined  as 

MAPE  =  lj2\Ak-Pk\/Ak.  (10) 

71  k= 1 

where  Ak  is  the  actual  (experimental)  value,  Pk  is  the  theoretically 
predicted  value  and  n  is  the  number  of  experiments.  MAPE  value 
is  4.0%  for  this  figure. 

In  Fig.  4,  the  variation  of  dimensionless  electrical  conductivity 
with  mean  temperature  is  given.  It  is  seen  that  electrical  conduc¬ 
tivity  considerably  increases  with  decreasing  temperature  and  it 
takes  nearly  4.5  times  higher  values  at  100  I<  than  the  value  at 
300  K.  This  is  because  of  the  lower  probability  of  collision  of  elec¬ 
trons  and  holes  with  lattice  atoms  and  impurities  at  lower  temper¬ 
atures.  MAPE  value  of  the  solid  line  is  4.5%  for  this  figure. 

The  measured  values  of  Seebeck  coefficient  and  electrical  con¬ 
ductivity  of  a  single  Bi2Te3  p-n  couple  at  300  K  are  otpn  =  345  pV/ 
°C  and  opn  =  28103  Q  m  respectively  while  they  are  given  as  otpn  = 
382  pV/°C  and  opn  =  30239  Q  m  in  literature  [19].  Therefore,  the 
relative  differences  between  the  measured  quantities  and  the  val- 


Fig.  3.  Variation  of  dimensionless  Seebeck  coefficient  with  mean  temperature  for 
AT  =  20  I<. 


Mean  temperature,  T(K) 

Fig.  4.  Variation  of  dimensionless  electrical  conductivity  with  mean  temperature 
for  AT  =  20  I<. 


ues  given  in  literature  are  less  than  10%  for  Seebeck  coefficient  and 
8%  for  electrical  conductivity. 


5.  Examination  of  maximum  power 


The  variations  of  the  material  properties  at  low  temperatures 
cause  the  change  of  behaviors  of  a  commercial  TEG  module.  Max¬ 
imum  power  output  can  be  calculated  by  using  the  experimentally 
measured  material  properties.  The  calculations  are  made  for  the 
same  temperature  difference  of  200  K,  which  is  the  common  oper¬ 
ating  temperature  difference  for  commercial  TEGs,  and  only  the 
mean  temperature  is  changed  between  100  K  and  375  K.  The 
dimensionless  form  of  maximum  power  output  is  defined  by  divid¬ 
ing  the  values  to  the  value  at  mean  temperature  of  300  K.  Maxi¬ 
mum  power  output  is  expressed  as 


Pmax(T,  AT) 


vie  [^Z^Tdrf 

4 Rj  a  Npn_  <5pn  rT+AT/2  __1 _Arr 

^  A  AT  Jf-AT/2  (Tpn(T)U1 


(11) 


Therefore  the  dimensionless  maximum  power  output  is  given 
by 


Pmax(T,  200) 


Pmax(T,  200) 
Pmax(300,  200)  • 


(12) 


It  should  be  noted  that  the  temperature  difference  is  kept  con¬ 
stant  as  AT  =  200  K  while  the  mean  temperature  is  changed.  Max¬ 
imum  power  output  of  TEG  decreases  with  decreasing  mean 
temperature  and  there  is  a  sharp  descent  of  power  output  after 
the  mean  temperature  of  200  K  as  it  seen  in  Fig.  5.  To  generate 
the  same  electrical  power  output  at  mean  temperature  of  211  K, 
which  corresponds  to  LNG  evaporation  application  (TL  =  -162  °C, 
Th  =38°C),  it  is  necessary  to  use  1.2  times  more  TEG  modules 
working  at  mean  temperature  of  41 1  K,  which  corresponds  to  stan¬ 
dard  application  (TL  =  38  °C,  TH  =  238  °C)  for  Bi2Te3  based  TEGs. 
Similarly,  for  evaporation  of  liquid  nitrogen  application  (TL  = 
-196  °C,  Th=4  °C)  it  is  necessary  to  use  1.45  times  more  TEG  mod¬ 
ules  in  comparison  with  their  standard  applications.  Furthermore, 
there  is  a  very  smooth  variation  of  dimensionless  power  output  be¬ 
tween  250  K  and  350  K.  It  is  nearly  equal  to  unity  for  this  temper¬ 
ature  interval.  Therefore  250  K  seems  a  critical  mean  operating 
temperature  for  this  particular  Bi2Te3  module  since  there  is  a  sharp 
descent  below  this  point. 
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Mean  temperature,  t(K) 

Fig.  5.  Dimensionless  maximum  power  output  versus  mean  temperature  for 
AT  =  200  K. 


6.  Conclusion 

TEGs  based  on  Bi2Te3  semiconductors  are  commercially  avail¬ 
able  and  the  cheapest  versions  of  thermoelectric  generators.  Their 
typical  low  and  high  operating  temperatures  are  50  °C  and  250  °C 
respectively.  In  this  work,  temperature  dependencies  of  Seebeck 
coefficient  and  electrical  conductivity  are  experimentally  deter¬ 
mined  for  low  temperature  range  100-375  K.  Some  expressions 
based  on  the  experimental  measurements  are  given  for  the  tem¬ 
perature  dependency  of  material  properties.  It  is  observed  that 
all  these  results  are  in  a  good  agreement  with  the  literature 
[14,16]  and  [19].  Furthermore,  by  using  the  experimentally  mea¬ 
sured  material  properties,  maximum  electrical  power  output  of  a 
Bi2Te3  based  TEG  module  is  analyzed  for  low  temperature  region. 
It  seems  that  mean  operating  temperature  of  250  K  constitutes  a 
critical  temperature  for  this  particular  thermoelectric  generator 
module.  Below  this  point,  it  is  necessary  to  use  more  TEG  modules 
for  the  same  power  output.  The  number  of  excess  thermoelectric 
modules  rapidly  increases  especially  below  200  K.  The  results  can 
be  used  for  the  low  temperature  applications  of  Bi2Te3  based  TEGs. 
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